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SUMMARY/OVERVIEW 


The  objectives  of  the  program  were  two  fold,  namely  (1)  to  identify  the  controlling  chemistry  of 
aeropropulsion  fuels  and  develop  detailed  and  reduced  mechanisms  describing  their  oxidative 
reaction  pathways  and  rates,  and  (2)  to  understand  and  quantify  the  unstable  modes  of 
combustion  and  their  influence  on  flame  extinction.  Concerning  Objective  1,  laminar  burning 
velocities  of  C,  to  C8  alkanes  with  air,  and  blends  of  butane,  ethylene  and  benzene  with  air,  were 
experimentally  measured.  These  results  were  used,  together  with  other  chemical  kinetic 
information  towards  the  development  of  a  comprehensive  C3-mechanism.  A  mathematical 
theory  based  on  computational  singular  perturbation  was  formulated,  allowing  for  the  systematic 
reduction  of  detailed  reaction  mechanism  to  simpler  ones  without  loss  of  the  comprehensiveness 
of  the  kinetic  description.  Concerning  Objective  2,  the  steady  and  pulsating  propagation  of 
flames  were  analytically  and  computationally  studied,  at  atmospheric  and  elevated  pressures. 
Results  showed  that  the  onset  of  pulsation  facilitates  flame  extinction,  that  positive  and  negative 
stretch  respectively  promote  and  retard  its  onset,  that  radiative  loss  promotes  flame  front 
instability  and  extinction,  and  that  there  is  strong  coupling  between  flame  front  stability  and  the 
intrinsic  chemical  chain  mechanisms.  These  results  are  expected  to  be  useful  to  the  general 
interests  of  AFOSR  in  the  fundamental  and  practical  issues  of  flame  dynamics  and  chemical 
kinetics,  turbulent  combustion,  radiative  hear  transfer,  flame  extinction,  stabilization, 
flammability,  and  supersonic  combustion. 


2 


ACCOMPLISHMENT 


Highlights  of  the  accomplishments  can  be  found  in  the  annual  reports  submitted  to  the  program 
director,  as  well  as  the  journal  and  conference  papers  that  have  appeared  in  print.  Thus  only  a 
brief  summary  of  these  works  is  mentioned  in  the  following. 

Development  of  Comprehensive  Kinetic  Mechanisms  for  Hydrocarbon  Oxidation 

Recognizing  the  lack  of  a  comprehensive  chemical  kinetic  mechanism  for  aeropropulsion  fuels  in 
general,  and  the  significantly  reduced  residence  time  and  hence  heightened  performance 
sensitivity  to  chemical  kinetics  in  supersonic  combustion  in  particular,  an  experimental- 
computational  effort  was  initiated  for  the  development  of  such  a  mechanism  for  hydrocarbon 
fuels  up  to  C8.  The  basic  requirements  for  this  mechanism  include:  (a)  fundamental  input  of 
thermochemistry  and  elementary  reaction  kinetics;  and  (b)  predictive  capability  of  the  wide 
spectrum  of  combustion  phenomena,  from  the  response  of  homogeneous  systems  such  as  shock 
tube  and  auto-ignition  delays,  to  the  ignition,  steady  burning,  and  extinction  of  premixed  and 
diffusion  flames  for  which  diffusive  transport  is  present,  and  ultimately  to  the  characteristics  of 
pollutant  formation.  The  reaction  mechanism  should  also  be  applicable  to  a  wide  variety  of 
hydrocarbon  fuels  and  should  incorporate  any  pressure  effects  on  the  combustion  responses.  To 
achieve  this  goal,  experiments  were  conducted  to  acquire  fundamental  combustion  data  such  as 
the  laminar  burning  velocities  and  the  extinction  of  stretched  premixed  and  diffusion  flames, 
while  ab  initio  quantum  mechanical  calculations  were  also  performed  to  determine  the  rate 
coefficients.  To  bring  the  fundamental  understanding  of  reaction  kinetics  to  a  predictive  model, 
systematic  optimization  of  the  reaction  rate  parameters  within  their  uncertainties  was  performed, 
and  the  optimized  model  validated  against  an  even  wider  range  of  combustion  data.  Lastly, 
attempts  were  initiated  to  reduce  the  detailed  reaction  mechanisms  to  smaller  mechanisms  that 
are  still  quantitatively  predictive  but  are  sufficiently  simple  for  implementation  in  large-scale 
modeling. 

During  this  period  of  funding,  significant  progress  was  made  on  the  combustion  kinetics 
of  low  molecular-weight  fuels,  including  acetylene,  ethylene,  propyne,  allene,  propene,  propane, 
and  butadiene.  The  reaction  kinetics  of  these  fuels  constitutes  a  critical  subset  of  a 
comprehensive  mechanism  of  aeropropulsion  fuels.  In  addition  to  the  study  of  small 
hydrocarbons  and  recognizing  the  lack  of  reliable,  fundamental  combustion  data  for  alkane  and 
alkene  fuels  in  general,  a  comprehensive  experimental  study  on  the  laminar  burning  velocities  of 
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Cj  to  C8  hydrocarbons,  as  well  as  some  fuel  blends,  was  also  conducted.  These  data  are  critical 
for  the  present  and  future  tasks  of  mechanism  development. 

Theoretically,  high-level  ab  initio  calculations  were  performed  for  a  number  of  reaction 
systems  pertinent  to  C2  and  C3  combustion.  The  reaction  rate  constants  were  determined  by  the 
Rice-Ramsperger-Kassel-Marcus  (RRKM)  theory.  Numerically,  sensitivity  analysis  was 
conducted  using  detailed  reaction  mechanisms,  which  were  validated  against  a  variety  of 
experimental  data,  including  fuel  pyrolysis  and  oxidation  in  shock  tubes,  laminar  burning 
velocities,  and  detailed  species  time  and  spatial  profiles  in  static  reactors  or  in  combusting  flows. 

Progress  made  in  this  endeavor  is  reported  in  the  following. 

1.  Combustion  Chemistry  of  Propane:  A  Case  Study  of  Reaction  Mechanism  Optimization 

Detailed  reaction  mechanisms  describing  hydrocarbon  combustion  have  hierarchical 
structures  with  H,  and  CO  chemistry  at  the  base,  supplemented  by  necessary  reaction  channels  of 
larger  chemical  species.  While  this  strategy  provides  rational  organization  of  combustion 
chemistry,  it  is  unknown  whether  parameter  sets  derived  by  optimizing  smaller-hydrocarbon 
chemistry  are  secure  foundations  for  extension  to  larger-hydrocarbon  oxidation  models.  This 
question  was  addressed  in  the  present  work  by  combining  the  newly  optimized  GRI-Mech  3.0  for 
methane  and  ethane  combustion  with  the  C3H  chemistry  developed  by  us  previously.  Model 
predictions  were  compared  against  methane,  ethylene,  ethane,  propyne,  allene,  propene,  propane 
laminar  burning  velocities  and  shock  tube  ignition  delay  data.  It  was  found  from  extensive 
sensitivity  testing  that  the  coupling  between  the  C.H^  chemistry  and  the  C,-C2  kinetic  subset  was 
much  stronger  than  anticipated.  The  Cl  combustion  data  could  not  be  satisfactorily  modeled 
using  the  C,-C2  kinetic  subset  without  re-optimizing  some  of  the  C<3  kinetic  parameters.  These 
parameters  have  stronger  influences  on  C3Hx  combustion  than  the  C3  parameters.  Parameter 
optimization  analogous  to  the  GRI-Mech  3.0  optimization  was  carried  out  using  essentially  the 
same  methods  as  were  used  in  the  development  of  C,-C2  kinetic  subset,  including  21  optimization 
targets,  of  which  9  were  ignition  delays  and  12  were  atmospheric  pressure  laminar  burning 
velocities.  A  reasonable  match  to  the  whole  target  set  could  be  obtained,  without  deteriorating 
the  agreement  between  experiment  and  calculation  for  the  methane  and  ethane  combustion  data, 
by  re-optimizing  9  parameters  from  the  C,-C2  kinetic  subset. 

The  implication  of  this  study  is  significant  in  that  any  reaction  mechanisms  proposed 
previously  and  many  more  to  come  represent  just  individual  points  within  the  uncertainty  space 
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of  a  given  kinetic  set.  Maximizing  the  likelihood  that  a  kinetic  model  can  predict  a  combustion 
property  beyond  the  regime  within  which  the  model  is  validated  can  only  be  achieved  by 
comprehensive  optimization.  In  this  way  the  size  of  the  uncertainty  space  is  reduced  only  to  an 
extent  as  being  inconsequential  to  the  combustion  of  a  given  set  of  target  fuels.  The  second 
implication  of  the  study  is  that  kinetic  model  optimization  in  a  multi-parameter,  -target,  and  -fuel 
space  is  entirely  possible  as  will  be  pursued  in  our  research  program.  Through  this  study,  we  also 
developed  a  tentatively  optimized  combustion  model,  which  predicts  very  well  a  wide  variety  of 
C,-C3  fuel  combustion. 

2.  Propene  Pvrolvsis  and  Oxidation  Kinetics  in  Flow  Reactor  and  Laminar  Flames 

The  pyrolysis  and  oxidation  of  propene  were  studied  experimentally  in  an  atmospheric 
flow  reactor.  Species  profiles  were  obtained  in  the  intermediate  to  high  temperature  range 
(-1200  K)  for  lean,  stoichiometric,  rich,  and  pyrolytic  conditions.  Laminar  burning  velocities  of 
propene/air  mixtures  were  also  determined  over  an  extensive  range  of  equivalence  ratios,  at  room 
temperature  and  atmospheric  pressure,  using  the  counterflow  twin  flame  configuration.  A 
detailed  chemical  kinetic  model  consisting  of  465  reactions  and  71  species  was  used  to  describe 
the  high  temperature  kinetics  of  propene,  propyne,  allene,  and  propane.  It  was  shown  that  the 
kinetic  model  could  accurately  predict  a  wide  range  of  combustion  data  for  these  fuels,  including 
laminar  burning  velocities,  speciation  in  flow  reactors,  and  ignition  in  shock  tubes.  Notable 
uncertainties  in  the  reaction  kinetics  of  these  fuels  are  identified  and  discussed. 

3.  RRKM.  Flow  Reactor,  and  Detailed  Kinetic  Modeling  Study  of  Propyne  Pyrolysis  and 
Oxidation 

The  oxidation  of  propyne  was  studied  experimentally  in  an  atmospheric  flow  reactor  and 
in  laminar  premixed  flames.  Species  profiles  were  obtained  for  propyne  oxidation  experiments 
conducted  in  the  Princeton  Turbulent  Flow  Reactor  in  the  intermediate  to  high  temperature  range 
(-1170  K)  for  lean,  stoichiometric,  and  rich  conditions.  A  detailed  chemical  kinetic  model  of 
high-temperature  propyne  oxidation,  consisting  of  437  reactions  and  69  species,  was  developed. 
It  was  shown  that  this  kinetic  model  predicts  reasonably  well  the  flow-reactor  and  flame-speed 
data  determined  in  the  present  study,  and  the  shock-tube  ignition  data  available  in  the  literature. 
The  remaining  uncertainties  in  the  reaction  kinetics  of  propyne  oxidation  were  discussed  and  an 
important  part  of  the  reaction  kinetics  was  examined. 
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The  pressure-dependent  rate  coefficients  for  several  reactions  relevant  to  propyne 
pyrolysis  were  determined  with  ab  initio  quantum  mechanical  calculations  and  Rice- 
Ramsperger-Kassel-Marcus  (RRKM)  analyses.  These  include  the  mutual  isomerization  of 
propyne  and  allene,  the  chemically  activated  reactions  of  propyne  and  allene  with  the  H  atom  and 
of  acetylene  with  methyl  on  the  C3H5  potential  energy  surfaces.  Propyne  pyrolysis  was  also 
experimentally  studied  in  a  flow  reactor  at  1210  K  and  1  atm.  A  detailed  reaction  mechanism, 
employing  the  current  RRKM  rate  coefficients,  was  shown  to  accurately  predict  the  experimental 
acetylene  and  methane  profiles  determined  in  the  flow  reactor  and  literature  shock-tube  data  of 
propyne  and  allene  pyrolysis  up  to  1500  K. 

4.  Quantum  Mechanical  and  Kinetic  Modeling  Studies  of  the  Initiation  Reactions  in  Acetylene 
Oxidation  in  Shock  Tube 

The  reaction  kinetics  of  acetylene  underlies  all  hydrocarbon  combustion  kinetics. 
Acetylene  oxidation  has  been  extensively  studied  in  shock  tubes.  A  large  amount  of  data  exists, 
which  can  be  used  to  verify  the  acetylene  subset  of  the  detailed  reaction  mechanism  for  aero- 
propulsion  fuels.  However,  validation  of  this  mechanism  subset  has  been  hampered  by  a  long¬ 
standing  and  unresolved  issue  in  acetylene  oxidation,  namely  the  nature  of  the  oxidation 
initiation  processes.  For  this  reason,  the  reaction  between  acetylene  and  molecular  oxygen  was 
analyzed  using  quantum  mechanical  calculations  and  kinetic  modeling  of  acetylene  oxidation  in 
shock  tubes.  Calculations  at  the  G2(B3LYP)  level  of  theory  showed  that  the  direct  attack  of 
molecular  oxygen  on  the  p  bond  in  acetylene  has  a  larger  energy  barrier  than  acetylene- 
vinylidene  isomerization,  such  that  this  isomerization  followed  by  the  reaction  of  vinylidene  with 
molecular  oxygen  is  the  energetically  favorable  initiation  reaction  of  acetylene  oxidation.  It  is 
further  shown  that  detailed  kinetic  models  of  acetylene  oxidation  including  this  initiation  process 
predict  well  the  experimental  shock  tube  ignition  delay  data.  Through  this  study,  it  is  the  first 
time  that  a  carbene  species  is  proposed  as  the  major  player  in  the  initiation  process  of  unsaturated 
hydrocarbon  oxidation.  It  remains  to  be  seen  whether  such  a  mechanism  is  generally  applicable 
to  the  oxidation  of  other  unsaturated  hydrocarbons  as  well. 

5.  Experiments  and  Detailed  Kinetic  Modeling  of  1.3-Butadiene  Oxidation  at  High  Temperatures 

The  oxidation  kinetics  of  1,3-butadiene  (1,3-C4H6)  is  of  considerable  importance  to  the 
hierarchical  development  of  the  kinetic  mechanisms  of  long-chain  aliphatics  and  aromatics 
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combustion.  Previously  a  comprehensive  kinetic  model  of  acetylene,  ethylene,  allene,  propyne, 
propene,  and  propane  combustion  was  developed.  This  model  can  predict  a  wide  variety  of 
combustion  responses  from  ignition  and  laminar  flame  propagation  to  the  overall  and  detailed 
structures  of  fuel  oxidation  and  combustion.  This  previous  model  was  extended  to  include  the 
kinetics  of  1,3-butadiene  pyrolysis  and  oxidation. 

Specifically,  the  high-temperature  kinetics  of  1,3-butadiene  oxidation  was  examined  with 
detailed  kinetic  modeling.  To  facilitate  model  validation,  flow  reactor  experiments  were  carried 
out  for  1,3-butadiene  pyrolysis  and  oxidation  over  the  temperature  range  1035-1185  K  and  at 
atmospheric  pressure,  extending  similar  experiments  found  in  the  literature  to  a  wider  range  of 
equivalence  ratios  and  temperatures.  The  kinetic  model  was  compiled  on  the  basis  of  an 
extensive  review  of  literature  data.  The  reaction  pathways  and  rate  constants  of  several  key 
reactions  were  examined  closely  by  a  detailed  analysis  of  previous  experimental  results,  by 
thermochemical  considerations,  and  by  Rice-Ramsperger-Kassel-Marcus  (RRKM)  calculations. 
Notably  these  reactions  include  the  attack  of  the  H  and  O  atoms  on  1,3-butadiene.  The  accuracy 
of  the  rate  constants  associated  with  these  reactions  proves  to  be  critical  to  the  overall  kinetics  of 
1,3-butadiene  pyrolysis  and  oxidation. 

The  kinetic  model  was  critically  validated  against  a  total  of  15  sets  of  experimental  data 
including  species  profiles  during  1,3-butadiene  pyrolysis  and  oxidation  in  flow  reactors  and 
shock  tubes,  and  1,3-butadiene  flames.  It  was  shown  that  the  kinetic  model  compiled  in  this 
study  is  capable  to  closely  predict  a  wide  range  of  high-temperature  oxidation  and  combustion 
responses. 

The  analysis  of  individual  reactions  and  the  reaction  pathway  based  on  the  current  kinetic 
model  revealed  a  viewpoint  of  1,3-butadiene  kinetics  significantly  different  from  previous 
studies.  Three  parallel  pathways  were  identified  to  be  relevant  to  1,3-butadiene  oxidation,  with 
the  chemically  activated  reaction  of  H  and  1,3-butadiene  to  produce  ethylene  and  the  vinyl 
radical  being  the  most  important  channel  over  all  experimental  conditions.  Extensive  sensitivity 
analyses  were  performed;  and  the  remaining  uncertainty  of  the  butadiene  chemistry  was 
discussed. 

6.  Determination  of  and  Fuel  Structure  Effects  on  Laminar  Burning  Velocities  of  Ch  to  C« 
Hydrocarbons 
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Laminar  burning  velocities  determined  by  using  the  counterflow  twin  flame  configuration 
were  compared  for  various  C,  to  C8  hydrocarbons,  including  alkanes,  alkenes,  alkynes, 
aromatics,  and  alcohols.  The  data  were  compared  over  an  extensive  range  of  equivalence  ratios 
at  room  temperature  and  atmospheric  pressure.  The  comparison  showed  that  the  laminar  burning 
velocities  of  normal  alkanes  are  close  throughout  the  entire  range  of  equivalence  ratios  studied, 
except  for  methane  whose  flame  speeds  are  consistently  lower.  The  more  unsaturated  the 
molecule  the  higher  the  burning  velocity  for  fuels  having  the  same  carbon  number  in  the  order  of 
alkanes  <  alkenes  <  alkynes.  Methyl  substitution  for  hydrogen  or  branching  reduces  the  flame 
speeds  for  both  alkanes  and  alkenes.  The  burning  velocities  of  large  saturated  cyclic  species 
(cyclohexane  and  cyclopentane)  are  close  to  those  of  their  normal  alkane  analogs. 

7.  Laminar  burning  Velocities  of  Ethvlene/n-Butane/Air  Mixtures 

The  objectives  of  the  present  study  were:  (a)  to  extend  previous  study  involving  single 
fuels  to  fuel  blends,  and  (b)  to  experimentally  determine  the  burning  velocities  by  using  Digital 
Particle  Image  Velocimetry  (DPIV). 

Regarding  the  first  objective,  it  is  noted  that  since  nearly  all  practical  fuels  are  multi- 
component,  it  is  important  to  develop  theoretical  and  semi-empirical  capabilities  to  predict  the 
flame  speeds  of  fuel  mixtures.  Furthermore,  existing  reaction  models  of  hydrocarbon 
combustion  were  mostly  developed  based  on  results  from  single  fuels  and  as  such  have  not  been 
adequately  verified  for  fuel  mixtures.  Finally,  it  is  also  of  interest  to  develop  semi-empirical 
mixing  rules. 

Mixtures  of  ethylene  and  n-butane  were  selected  for  study  because  ethylene  is  an 
important  intermediate  in  hydrocarbon  oxidation  while  n-butane  is  a  representative  n-alkane. 
The  experimental  data  of  the  laminar  burning  velocities  of  ethylene/air  and  n-butane/air  mixtures 
were  found  to  compare  well  with  the  calculated  values  obtained  by  using  the  kinetics  model 
developed  under  the  present  program.  Using  this  kinetics  model,  further  calculations  were 
performed  for  mixtures  of  ethylene  and  n-butane  in  the  ratios  of  2:1,  1:1,  and  1:2.  As  is 
reasonable  to  expect,  the  burning  velocities  of  the  fuel  mixtures  are  bounded  by  those  of 
ethylene/air  and  n-butane/air.  However,  results  for  the  mixtures  seem  to  be  weighted  towards  the 
n-butane  values.  The  same  biasing  is  also  exhibited  for  the  calculated  adiabatic  flame 
temperatures,  indicating  that  the  flame  temperature  instead  of  kinetic  coupling  could  have  the 
dominant  influence  on  the  burning  velocity.  The  biasing  could  simply  be  a  consequence  of  how 
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the  composition  of  mixture  is  defined.  The  present  mixing  ratio  is  the  conventional  one,  based 
on  the  molar  ratios.  There  are,  however,  other  definitions  that  could  be  more  relevant  physically. 
This  aspect  is  being  studied. 

Studies  on  Chemistry-Related  Problems 

8.  Effect  of  Inelastic  Collision  on  Binary  Diffusion  Coefficients  of  Free  Radicals 

Diffusion  is  a  key  factor  in  chemically  reacting  flows  in  which  the  propagation  rate  of  a 
flame  is  determined  by  diffusion  as  well  as  chemical  reaction  and  heat  release.  It  is  known  that 
the  sensitivity  of  flame  propagation  to  the  H-atom  diffusion  coefficient  is  as  large  as  that  of  the 
free-radical  chain  branching  reaction  H  +  02  — >  O  +  OH.  Despite  the  recent  progress  in  the 
description  of  H-atom  diffusion  coefficient,  the  molecular  theory  of  H-atom  diffusion  remains  to 
be  that  of  Chapman-Enskog,  in  which  elastic  collision  is  assumed.  On  the  other  hand,  it  is 
known  that  collision  of  the  H  atom  with  notable  combustion  species  including  H,,  CO,  C02,  02, 
C2H>,  and  C2H4  is  inelastic  at  combustion  temperatures.  At  very  high  temperatures,  say,  2500  K, 
nearly  a  quarter  of  collision  encounter  between  the  H  atom  and  N2  is  inelastic.  This  inelastic 
collision  can  be  best  described  by  the  formation  of  a  transient  chemical  bond,  whose  binding 
energy  is  of  the  order  of  103  times  larger  than  the  well  depth  of  a  typical  Lennard-Jones  potential. 
In  this  work,  we  examined  the  influence  of  inelastic  collision  on  diffusion  coefficients.  Our  aim 
was  to  obtain  an  order-of-magnitude  estimate  for  this  influence,  and  compare  it  to  the  influence 
of  reaction  rate  parameters.  In  this  way  the  uncertainty  resulting  from  the  diffusion  coefficients 
could  be  assessed. 

Diffusion  coefficient  was  examined  using  molecular  dynamics  simulation  and  the  Green- 
Kubo  formula,  with  pair  interaction  potential  energies  corresponding  to  collision  with  and 
without  the  elastic  assumption.  It  was  found  that  inelastic  collision  significantly  increases  the 
diffusion  coefficients  of  free  radical  species  at  temperatures  relevant  to  fuel  combustion. 
Inelastic  collision  with  a  potential  energy  curve  representative  of  molecule-radical  interaction 
may  cause  the  diffusion  coefficient  to  increase  by  as  much  as  50%  at  1500  K.  This  finding 
suggests  that  a  molecular  theory  beyond  the  Chapman-Enskog  equation  is  needed  to  describe 
accurately  free-radical  diffusion  in  laminar  flames. 

9.  Theory  of  Complex  Computational  Singular  Perturbation 
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The  method  of  Computational  Singular  Perturbation  (CSP)  for  the  analysis  and  reduction 
of  detailed  chemical  mechanisms  was  extended  to  the  complex  eigensystem.  Results  showed 
that  the  time  scales  of  chemical  species  change  dramatically  and  non-monotonically,  and  that 
oscillatory  modes  appear  frequently  in  large  chemical  reaction  mechanisms.  The  present  method 
was  then  employed  to  generate  reduced  mechanisms  for  hydrogen/air  and  methane/air  oxidation. 

Using  the  hydrogen/air  system  as  an  example,  it  was  found  that  the  number  of  complex 
modes  increases  with  decreasing  residence  time  and  hence  increasing  chemical  sensitivity,  and 
that  the  longer  the  residence  time  the  more  the  number  of  steady-state  species  and  hence  the 
smaller  the  mechanism.  The  validity  of  these  reduced  mechanisms  was  evaluated  based  on  the 
responses  of  the  perfectly  stirred  reactors  and  the  one-dimensional  planar  propagating  premixed 
flames.  Comparison  between  the  reduced  and  detailed  chemistries  over  a  wide  range  of 
pressures  and  equivalence  ratios  showed  good  agreement  on  the  flame  speed,  temperature  and 
structure. 

10.  Structure  and  Sooting  Limits  of  Diffusion  Flames  from  1  to  5  Atmospheres 

Previous  studies  on  the  structure  and  soot  formation  of  atmospheric,  non-sooting  and 
near-sooting  counterflow  ethylene  diffusion  flames  showed  that  the  experimental  soot  volume 
fractions  and  laser  scattering  signals  can  be  satisfactorily  predicted  by  a  detailed  kinetic  model  of 
soot  formation.  By  further  extending  the  study  to  include  the  influence  of  pressure  and  fuel 
dilution  on  the  critical  sooting  strain  rates,  it  was  identified  that,  when  the  flame  temperature  is 
held  constant,  the  density-weighted  strain  rate  at  the  sooting  limit  can  be  correlated  with  the 
system  pressure,  p,  and  fuel  mole  fraction,  XF,  according  to  p0Kp  ~  pXF05,  where  p0  is  the  mass 
density  of  the  cold  oxidizer  gas  mixture,  and  Kp  the  local  strain  rate  at  the  sooting  limit.  Hence, 
a  plot  of  p0Kp  versus  pXF0'5  defines  the  sooting  boundary  in  that  the  region  above  the  limit  is 
non-sooting,  and  that  below  it  is  sooting.  It  was  further  shown  that  such  a  correlation  could  be 
understood  by  plotting  pXF° 5  versus  the  peak  acetylene  partial  pressure,  which  again  exhibits  a 
linear  dependence.  Hence  a  linear  relationship  exists  between  the  sooting  limit,  p0Kp,  and  the 
partial  pressure  of  acetylene,  implying  that  a  larger  acetylene  partial  pressure  requires  stronger 
straining  to  prevent  soot  formation. 

The  sooting  limit  correlation  identified  highlights  the  practical  significance  of  such  a 
correlation,  in  that  it  provides  an  unambiguous  means  to  predict  whether  a  flame  environment 
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with  given  strain  rate,  pressure,  and  fuel  mole  fraction  would  result  in  a  local  sooting  condition 
in  a  complex  flame  environment  such  as  the  turbulent  diffusion  flame.  In  addition,  the  simplicity 
of  the  correlation  indicates  its  fundamental  importance,  which  is  yet  to  be  fully  understood 
through  kinetic  modeling  studies. 

The  present  work  extended  the  study  to  methane  and  propane  flames.  The  temperature 
and  major  species  concentration  profiles  were  measured  with  spontaneous  Raman  scattering,  and 
computationally  simulated  with  detailed  kinetics  and  transport.  Good  agreement  was  found 
between  the  experimental  data  and  the  computational  simulation,  hence  providing  benchmark 
data  for  further  studies.  It  was  also  shown  that  the  previously  developed  global  and  local  sooting 
limit  correlations  are  again  applicable,  respectively  relating  the  density- weighted  strain  rate  at  the 
sooting  limit  with  the  global  parameters  of  the  system  pressure  and  fuel  mole  fraction  in  the  fuel 
stream,  and  with  the  local  flame  parameter  of  the  peak  acetylene  partial  pressure.  This  again 
demonstrated  the  importance  of  acetylene  as  a  crucial  intermediate  in  soot  formation  in  flames. 
In  addition,  the  local  correlations  for  the  propane  and  ethylene  flames  collapse  into  a  single 
relation.  The  fact  that  the  correlation  for  the  methane  flames  is  different  from  those  of  propane 
and  ethylene  is  not  unexpected  considering  the  significantly  different  kinetic  features  of  methane 
from  those  of  other  fuels,  particularly  in  terms  of  the  detailed  reaction  kinetics  of  PAH  and  soot 
formation  from  methane,  which  is  yet  to  be  fully  understood. 

Unsteady  Flame  Dynamics 

11.  Steady  and  Pulsating  Propagation  of  Rich  Hvdrogen/Air  Flames  at  Atmospheric  and 
Elevated  Pressures 

Previous  computational  simulation  of  atmospheric  one-dimensional  planar  hydrogen/air 
flames  with  detailed  chemistry  and  transport  showed  that  pulsating  instability  develops  as  the 
rich  flammability  limit  is  approached.  Specifically,  for  <(>=73,  the  maximum  flame  temperature 
and  hence  flame  propagation  remain  steady.  However,  the  propagation  becomes  oscillatory  at 
((>=7.4,  changes  to  period  doubling  at  ((>=7.6,  and  extinguishes  at  (J>=7.8.  Since  the  steady-state 
extinction  turning  point  occurs  at  ((>=10.22,  the  result  showed  that  oscillation  reduces  the 
flammability  limit  defined  for  the  steady  flame. 

The  present  study  extended  the  previous  investigation  to  include  effects  of  elevated 
ambient  pressure,  for  pressures  ranging  from  1  to  20  atmospheres,  with  emphasis  on  the 
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influence  of  chain-mechanisms  in  H2/02  oxidation.  For  steady  burning,  with  increasing  pressure, 
the  burning  rate  first  increases,  then  decreases,  and  then  increases  again,  thereby  resulting  in  the 
occurrence  of  a  maximum  burning  rate  and  a  minimum  burning  rate.  It  was  further  found  that 
the  states  of  maximum  and  minimum  burning  rates  largely  correspond  to  those  of  the  crossover 
temperature.  Sensitivity  analysis  then  demonstrated  that  the  first  regime  is  dominated  by  the 
H+02,  0+H2,  and  OH+H2  fast-branching  sequence,  the  second  regime  is  significantly  influenced 
by  the  H+02+M  terminating  reaction,  and  the  third  regime  is  dominated  by  the  slow-branching 
H02  and  H202  chemistry.  It  is  significant  to  recognize  that  the  response  of  the  burning  rate, 
which  is  expected  to  be  fundamentally  affected  by  diffusion,  follows  that  of  the  explosion 
chemistry  of  the  homogeneous  system. 

Further  calculations  were  performed  for  the  states  of  steady  extinction,  the  onset  of 
pulsation,  and  the  extinction  of  the  pulsating  flame.  It  was  again  observed  that  pulsating 
extinction  occurs  at  lower  concentrations  than  those  of  the  steady  flames,  indicating  a  narrowing 
of  the  flammability  limit  when  considering  flame  unsteadiness.  Furthermore,  the  states  for  the 
onset  of  pulsation  are  very  close  to  those  of  the  maximum  burning  rate,  which  implies  that  the 
onset  of  pulsation  is  also  controlled  by  the  crossover  temperature.  The  cause  for  this  very  close 
coincidence  has  yet  to  be  identified,  although  it  must  be  somehow  related  to  the  slowing  down  of 
the  reaction  rate  by  the  H+02+M  reaction,  beyond  the  maximum  burning  state,  and  thereby  the 
increase  in  the  effective  activation  energy.  This  in  turn  favors  the  Sivashinsky  criterion  Ze(Le- 
1)>11  for  the  onset  of  pulsating  instability,  where  Ze  is  the  Zel’dovich  number  and  Le  the  Lewis 
number. 

The  study  also  found  that  a  mixture  that  is  pulsatingly  extinguishable,  and  hence 
nonflammable,  can  assume  a  steady  burning  mode  at  higher  pressures.  This  is  caused  by  the 
reduction  in  the  effective  activation  energy  for  the  weakly  branching  reactions  that  are  dominant 
under  high  pressures. 

12.  Pulsation  and  Flammability  Limits  of  Lean  Heptane/ Air  Flames 

The  unsteady  propagation  of  lean  heptane/air  planar  flames  in  the  doubly  infinite  domain 
was  computationally  simulated  with  detailed  chemistry  and  transport,  and  with  and  without 
radiative  heat  loss.  For  the  adiabatic  situation,  thermal-diffusive  pulsating  instability,  promoted 
for  large  values  of  the  Lewis  and  Zeldovich  numbers,  was  observed  to  develop  when  the  mixture 
becomes  sufficiently  lean.  With  subsequent  progressive  reduction  in  the  equivalence  ratio,  0,  the 
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mode  of  pulsation  changes  from  that  of  monochromatic,  to  period  doubling,  and  to  bursting 
characterized  by  bursts  of  high  burning  intensity  separated  by  long  period  of  dormancy.  The 
flame  nevertheless  does  not  extinguish.  This  behavior  is  similar  to  that  of  the  rich  hydrogen/air 
flame  studied  previously.  However,  when  heat  loss  is  considered,  the  onset  of  pulsation  is 
facilitated,  although  the  regime  in  <j)  for  sustained  pulsating  propagation  is  extremely  narrow,  and 
extinction  occurs  readily  when  the  amplitude  of  the  flame  temperature  oscillation  becomes  large 
enough  to  extinguish  the  flame  during  the  negative  phase  of  excursion.  The  abruptness  of 
transition  to  extinction  is  quite  different  from  that  of  the  rich  hydrogen/air  flame.  The  pulsating 
extinction  occurs  at  a  larger  </>  (=0.5063)  than  steady  extinction  (^=0.4762),  again  implying  that 
the  flame  extinguishes  in  the  pulsating  instead  of  the  steadily  propagating  mode,  and  that  the 
flammability  limit  is  accordingly  narrowed.  It  was  further  demonstrated  that  the  states  for  the 
onset  of  adiabatic  and  non-adiabatic  pulsations  can  be  accurately  estimated  by  respectively  using 
the  criterion  of  Sivashinsky  and  of  Joulin  and  Clavin,  provided  that  the  Lewis  number  used  is  the 
global  one,  extracted  from  the  response  of  stretched  flames.  This  therefore  demonstrates  that  the 
Lewis  number  is  a  property  of  the  flame  instead  of  the  unbumed  mixture  alone.  It  is  also 
suggested  that  the  lean  flammability  limits  of  mixtures  of  large  hydrocarbons  and  air  can  be 
estimated  by  using  the  criterion  of  Joulin  and  Clavin. 

13.  Effects  of  Stretch  on  Flame  Pulsation 

Effects  of  stretch  on  the  pulsating  instability  of  premixed  flames  were  investigated  via  the 
negatively  stretched  inwardly  propagating  spherical  flame  (IPF)  and  the  positively  stretched 
counterflow  flame  (CFF).  Computational  and  asymptotic  analysis  results  showed  that  pulsating 
instability  is  suppressed  by  the  former  and  promoted  by  the  latter.  In  particular,  it  was  shown 
that  for  a  given  rich  hydrogen/air  mixture  whose  one  dimensional,  freely  propagating  flame  is 
pulsatingly  unstable,  the  IPF  initially  propagates  at  the  laminar  burning  velocity  when  the  flame 
radius  is  large.  Oscillation  subsequently  develops,  and  is  then  amplified,  damped,  and  eventually 
suppressed  as  the  flame  propagates  inward  and  the  magnitude  of  stretch  increases.  For  the 
positively  stretched  rich  hydrogen/air  CFF,  oscillation  is  initiated  at  an  equivalence  ratio  much 
smaller  than  the  one-dimensional  rich  threshold.  Furthermore,  the  critical  strain  rate  leading  to 
pulsation  is  smaller  than  the  corresponding  static  extinction  limit,  implying  that  the  flame 
extinguishes  in  the  pulsating  instead  of  the  steadily  propagating  mode  such  that  the  flammable 
range  is  accordingly  narrowed.  In  addition,  it  is  seen  that  the  pulsating  flames  are  quasi-steady  in 
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nature  in  that  the  period  of  oscillation  is  larger  than  the  characteristic  flame  time.  As  such,  the 
unsteady  flame  cannot  recover  once  the  instantaneous  flame  temperature  is  reduced  below  the 
corresponding  steady-state  extinction  temperature. 

The  computationally  determined  pulsation  limits  using  realistic  chemistry  and  transport 
were  also  found  to  agree  well  with  the  asymptotic  results  based  on  simplified  chemistry  and 
transport,  provided  that  the  global  activation  energies  and  Lewis  numbers  are  appropriately 
extracted  from  laminar  flame  responses. 

14.  Oscillatory  Extinction  of  Diffusion  Flames 

The  transient  behavior  of  spherical  diffusion  flames  was  computationally  studied  in  both 
the  high-velocity  and  radiation-induced  limit  of  the  isola  response  of  flame  extinction. 
Oscillatory  instability  was  observed  near  both  steady-state  extinction  limits,  with  the  oscillation 
amplitude  growing  until  it  becomes  large  enough  to  extinguish  the  flame.  For  the  hydrogen  and 
methane  flames  studied,  oscillation  always  leads  to  extinction,  although  the  characteristics  of 
oscillation  are  qualitatively  different.  Specifically,  the  methane  flame  exhibits  large  amplitude 
temperature  fluctuations  at  a  frequency  of  about  0.35  Hz,  while  the  hydrogen  flame  has  small 
amplitude  fluctuation  with  a  frequency  of  about  60  Hz.  The  oscillations  in  the  methane  flame 
appear  to  be  characteristic  of  those  studied  based  on  the  diffusion  flame  structure,  while  those  of 
hydrogen  flames  are  reminiscent  of  those  observed  in  rich  hydrogen/air  premixed  flames. 

15.  Radiation  Induced  Flamefront  Instability 

Unsteady  propagation  and  instability  of  radiative  stretched  premixed  flames  at  Lewis 
numbers  close  to  unity  were  numerically  investigated  using  helium-diluted  counterflow  methane 
flames.  The  objective  was  to  understand  the  stability  of  the  three  near-limit  and  sub-limit 
multiple  flame  regimes,  namely  the  near  stagnation  flame,  the  weak  flame,  and  the  distant  flame. 
It  was  shown  that  both  the  near-limit  and  sub-limit  near-stagnation  flames  are  stable  close  to  the 
flammability  limit,  and  that  a  new  mode  of  oscillation  instability  exists  at  equivalence  ratios 
sufficiently  larger  than  that  at  the  flammability  limit.  In  addition,  stability  of  the  weak  and  the 
distant  flames  was  demonstrated,  while  unsteady  mutual  transitions  between  the  near-stagnation 
flame  and  the  distant  flame,  and  the  jump  from  the  weak  flame  to  the  normal  flame,  were  found. 
A  comparison  of  the  present  prediction  with  the  microgravity  experimental  observation  showed 
qualitative  agreement  on  the  oscillation  frequency  and  the  regime  for  the  occurrence  of 
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oscillation.  Calculations  for  Lewis  numbers  larger  and  smaller  than  unity  showed  that  the 
radiation  induced  instability  is  a  physically  intrinsic  phenomenon  for  radiative  stretched 
premixed  flame. 
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